molecular dynamics investigation of dislocation pinning by a nanovoid in copper 
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Interactions between an edge dislocation and a void in copper are investigated by means of molec- 
ular dynamics simulation. The depinning stresses of the leading partial and of the trailing partial 
show qualitatively different behaviors. The depinning stress of the trailing partial increases logarith- 
mically with the void radius, while that of the leading partial shows a crossover at 1 nm above which 
two partials are simultaneously trapped by the void. The pinning angle, which characterizes the 
obstacle strength, approaches zero when the void radius exceeds 3 nm. No temperature dependence 
is found in the critical stress and the critical angle. This is attributed to an absence of climb motion. 
The distance between the void center and a glide plane asymmetrically affects the pinning strength. 

PACS numbers: 61.80.Az, 62.20.Fe, 61.72.Qq 



I 

O 

o 



> 

ON 

O 

-I— > 



I 

O 

o 



X 



I. INTRODUCTION 

Voids are ubiquitous in irradiated metals and act as 
obstacles to dislocation motion as well as other radiation- 
induced defects: e.g. stacking fault tetrahedra or helium 
bubbles. Those obstacles results in the increase of yield 
stress and therefore play an important role in irradia- 
tion hardening. To investigate the extent of hardening 
by those obstacles, there has been a model in which a 
dislocation is regarded as a continuous line with a con- 
stant line tension. This is referred to as the uniform line 
tension model. In the presence of obstacles, a dislocation 
is fixed to form a cusp at an obstacle. The pinning angle 
4> is defined as the angle between two tangent vectors at 
a cusp. (See FIG.nj. Then the restoring force to make a 
dislocation straight is written as 27Cos(0/2), where 7 de- 
notes the line tension. We assume that a dislocation can 
penetrate an obstacle when the restoring force exceeds 
the critical value. Since 7 is a constant, this condition 
is equivalent to < 0c, which we call the critical angle. 
Note that stronger obstacles have smaller critical angles. 
A dislocation bows out to form an arc between two ob- 
stacles until the pinning angle reaches its critical value. 

For a periodic array of obstacles whose spacing is L, 
the critical resolved shear stress (CRSS) Tc above which 
a dislocation can penetrate the array of obstacles is rep- 
resented by 
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FIG. 1: A cusp formed at an obstacle. The angle 
two tangential vectors is called the pinning angle. 



where h denotes the Burgers vector length of the dislo- 
cation. (The line tension 7 is given by the elastic theory 
and is often written as G6^/2, where G represents the 
shear modulus.) 

In more realistic situations, obstacles are distributed 
randomly on glide planes and the randomness plays a 
crucial role in dislocation dynamics. In order to incorpo- 
rate this effect. Foreman and Makin performed computer 
simulation of the dislocation motion on a glide plane with 
randomly distributed obstacles which have the same crit- 
ical angle 0. They found that the dislocation propaga- 
tion has two qualitatively different modes depending on 
the critical angle. For obstacles of a small critical an- 
gle (i.e. strong obstacles), the dislocation propagation 
resembles dendritic growth while for a large critical 
angle (i.e. weak obstacles) the global form does not sig- 
nificantly deviate from the straight line. Also Tc are well 
described by 
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where L denotes the inverse of the square root of the 
areal density of obstacles on a glide plane 

In this context, the critical angle <j)c, which charac- 
terizes the obstacle strength, is an important parameter 
to discuss the extent of hardening. However, estimation 
of the critical angle is not an easy task, because it in- 
volves the core structure of dislocations. In this regard, 
extensive molecular dynamics (MD) simulations on the 
interaction between a dislocation and radiation-induced 
obstacles have been performed: e.g. stacking fault tetra- 
hedra Q, interstitial Frank loops voids and copper 
precipitates in bcc iron 0, 11,13 ■ 

On the other hand, there are some experimental at- 
tempts to determine the critical angle utilizing transmis- 
sion electron microscopy (TEM) They seem to 
be promising but still should be complemented by MD 
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FIG. 2: Scliematic of tlie system. Periodic boundary condi- 
tions are assigned to tlie x and tlie y directions. Tlie system 
contains approximately 0.55 million atoms. 

simulations in terms of the spatial resolution. For exam- 
ple, subnanometer obstacles cannot be seen by TEM. 

In this paper, along the line of the above computa- 
tional and experimental studies, we wish to estimate the 
critical angle and the critical stress regarding voids in 
copper. Especially, effects of the dissociation, temper- 
ature, and the distance from the void center to a glide 
plane are studied in detail. Note that we focus on an 
edge dislocation here. Results on a screw dislocation will 
be presented elsewhere. 

This paper is organized as follows. In the section^ we 
introduce the computational model; fee copper including 
an edge dislocation and a void. In the section IIIII the 
nature of the critical stress is investigated. Especially, 
void size dependence and the temperature dependence is 
discussed. The section Hvl deals with the critical angle, 
which is often measured by experiments (mainly TEIVI 
observation) to determine the obstacle strength. The dis- 
cussion enables us to compare the IMD simulation with 
the experiments. In the section Ivl effects of the distance 
between the void center and a glide plane are investi- 
gated. In the section IVII we discuss pinning strength 
of deformed voids in the context of dislocation channel- 
ing and plastic flow localization. The last section IVIII is 
devoted to discussions and conclusions. 



II. THE IVEODEL 

A. the geometry 

We treat fee copper in this paper. As for the inter- 
atomic potential, we adopt the embedded-atom method 
of Finnis-Sinclair type jl2| and choose the parameters 
according to Ackland et al. The lattice constant 

a = 3.615 A. 

The schematic of our system is shown in FIG. |2 The 
x, 2/, and z axes are taken as the [112], [110], and [111] 
directions, respectively. The length of each dimension is 
23nm, 23nm, and 15 nm. Periodic boundary conditions 
are employed in the x and the y directions. That is, we 
consider dislocations of infinite length in the x direction 



are periodically located in the y direction. 

Note that we have the surface only in the z direction, 
both for z > and for z < 0. Following three atomic 
layers of [111] next to the lowest surface (z < 0) are "the 
fixed layers" where velocities of the atoms always vanish. 
Similarly, three atomic layers of [111] next to the upper 
surface (z > 0) are "the moving layers" where velocities 
of the atoms are not given by the integration of force 
acting on them, but are given as a constant to cause 
the shear stress. Namely, the strain rate is the control 
parameter: not the shear stress. The moving surface 
moves to the —y direction: i.e. [110]. 

In order to introduce a void, atoms whose barycentric 
positions are in the spherical region {x±L/2)'^+y'^ + z^ < 

are removed, where r denotes the radius of the void. 
To introduce an edge dislocation, atoms that belong to 
one (110) half plane (z < 0) are removed and the rest of 
atoms are displaced by the strain field calculated from 
the elastic theory. This procedure produces a perfect 
edge dislocation whose Burgers vector is a/2[110]. How- 
ever, a perfect dislocation in an fee crystal is energeti- 
cally unstable to split into two partial dislocations whose 
Burgers vector length 6 is a /V6 = 1.48 A. 

^[110] ^^[211] + ^[121]. (3) 

Since we wish to prepare two partial dislocations and a 
void in an initial system, atoms are suitably shifted from 
the original position by the steepest descent method in 
order to realize the dissociation. In addition, since our 
system consists of the periodic array of dislocations due 
to the periodic boundary conditions in the x direction, 
this procedure also incorporates the excess strain field 
caused by the next dislocations. After certain time steps, 
the perfect dislocation dissociates to yield two partial 
dislocations separated by approximately 4 nm. 

Temperature is fixed to be 300 K in this paper, ex- 
cept for the section IIII CI where the temperature ef- 
fect on the CRSS is investigated. Velocities of atoms 
are given by random numbers which obey the IMaxwell- 
Boltzmann distribution. After a relatively short time re- 
quired for phonon relaxation, "the moving layer", which 
is explained above, begins to displace to cause strain. 

B. the strain rate 

The strain rate e is an important parameter in IVID sim- 
ulations of dislocations. In this paper, we set e = 8 x 10^ 
[sec~^]. Although it seems still unrealistically fast de- 
formation, the strain rate in MD simulation should not 
be directly compared with the macroscopic (or experi- 
mental) strain rate because the macroscopic strain rate 
involves only the average dislocation velocity. Namely, 
both spatial and temporal fluctuations in dislocation ve- 
locity are neglected. The complete correspondence of mi- 
croscopic dislocation velocity to the macroscopic strain 
rate is not clear at all unless we know the statistical prop- 
erty of space-time fluctuation in dislocation motion. 
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FIG. 3: Successive snapshots of a depinning process: (a) trap- 
ping of the leading partial, (b) just before the depinning of 
the leading partial, (c) just before depinning of the trailing 
partial, (d) depinning of the trailing partial. The void radius 
is 1 nm. To visualize the void and the dislocations, atoms that 
have 12 nearest neighbors are omitted. (Atoms which form 
dislocations have 11 or 13 nearest neighbors, and the number 
of nearest neighbors of void surface atoms is less than 12.) 



Calculation of the shear stress is noteworthy. We de- 
fine the shear stress as the forces m the y direction acting 
on the unit area of the moving layers and the fixed lay- 
ers. Hovifever, this microscopic definition of the shear 
stress shows a large thermal fiuctuation. In addition, 
the inertial motion of the dislocation due to the high 
strain rate may enable depinning at lower stress 0. In 
order to reduce these effects, simulation is performed 
twofold. Namely, the representative point in the phase 
space (spanned by the positions and the momenta of all 
atoms) are recorded every 4.6 ps. We then take the point 
where depinning just begins and restart the simulation 
with e = 0. In this relaxation process, the inertial ef- 
fect is ruled out and the thermal fiuctuation in the shear 
stress is averaged out. Taking this relaxation process into 
account, the average strain rate becomes approximately 
10^ [1/sec]. 



III. BEHAVIORS OF THE CRITICAL 
RESOLVED SHEAR STRESS 

A. temporal behavior 

First, we track the time evolution. Under the shear 
stress caused by the boundary condition, two partial dis- 
locations glide towards the [110] direction: the —y direc- 
tion. The Burgers vectors of the leading and the trailing 
partials are [211] and [121], respectively. Snapshots of a 
pinning-depinning process are shown in FIG.|21 

Note that there are two depinning processes of the 
leading and the trailing partials. In the stress-strain re- 
lation (FIG. El, we can see two peaks which correspond 
to the each depinning process. We remark that the crit- 
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FIG. 4: Stress-strain relation. Two arrows indicate the depin- 
ning points of two partials. Two circles which indicate sudden 
stress drops result from the attractive interaction between the 
dislocations and the void. 



ical stress for the leading partial is always larger than 
that of the trailing partial. This is because the leading 
partial has already escaped from the void and keeps glid- 
ing while the trailing partial is pinned by the void. Then 
the width of the stacking fault ribbon extends to yield at- 
tractive force between the partials. Namely, depinning of 
the trailing partial is assisted by the attractive force from 
the leading partial. Also, note that the movement of the 
leading partial leads to the stress relaxation, which ap- 
pears in the change of modulus in the stress-strain curve 
in FIG. El That is, the modulus is 30 GPa when the lead- 
ing partial is pinned (0.002 < e < 0.006), whereas the one 
after the leading partial is unpinned (0.008 < e < 0.012) 
is 13 GPa. 



B. void size dependence 

Then the critical stress is calculated for various voids 
of different radii (from 0.3 to 2.5 nm). Wc measure the 
depinning stresses both for the leading partial and for 
the trailing partial. An interesting feature arises from 
the comparison of the both stresses, which are shown in 
FIG. 13 We can see that they have different tendencies 
with respect to the void radius r. The depinning stress of 
the trailing partial shows the well-known logarithmic de- 
pendence, while that of the leading partial (i.e. the yield 
stress) shows a crossover around r ~ 1 nm. We will dis- 
cuss the difference more quantitatively in this subsection. 

The depinning stress of the trailing partial can be de- 
scribed by the following relation. 
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where T and B denote arbitrary constants. The best fit is 
realized by letting T = 37 MPa and S = 0.14 nm. Note 
that this logarithmic dependence has also been found in 
a continuous model by Scattergood and Bacon [l5j. and 
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FIG. 5: Void radius dependence of the depinning stress for 
the leading partial (x) and for the trailing partial (+)• The 
dashed line denotes Eq. ^ with T = 37 MPa and B = 0.14 
nm. The dotted line denotes Eq. (El, which is the result from 
a continuous model calculation 
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also in the context of the Orowan mechanism bmce 
their computation deals with only the critical stresses 
for pure edge dislocations and screw dislocations, it does 
not quantitatively apply to the partial dislocations. At 
least, it can describe the qualitative behavior regarding 
the trailing partial. 

On the other hand, the crossover of the yield stress 
cannot be explained in the above context. It seems to be 
describable by piecewise logarithmic behavior; 
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Note that the extent of hardening becomes much greater 
when the void radius exceeds 1 nm. Although we can- 
not show the definite explanation of this phenomenon, a 
plausible reason lies in the interaction between two par- 
tials. When the void radius is larger than 1 nm, the two 
partials are simultaneously trapped by the void. Then 
the depinning of the leading partial may be significantly 
affected by the trailing partial, since the distance between 
them is very close near the void (a few Burgers vectors) . 
Note that, for r = 0.3 nm and 0.5 nm, depinning of 
the leading partial takes place before the trapping of the 
trailing partial; i.e. the depinning of the leading partial 
is less affected by the trailing partial. We also remark 
that a similar crossover due to the interaction between 
dislocations is observed in the critical stress of dislocation 
nucleation on the void surface p^ . 

Regardless of the void radius, the depinning of the 
trailing partial takes place when the leading partial is 
completely escaped from the void. Although the trail- 
ing partial interacts with the leading partial via stacking 
fault ribbon and the elastic strain, it is not as strong as 
the one in the depinning process of the leading partial. 
Therefore, it is reasonable that the unpinning stress for 
the trailing partial can be described in the framework of 
|15| where a single dislocation involves. 




FIG. 6: Snapshots of dislocations (the trailing partials) just 
before depinning. The trailing partial dislocation at 100 K 
is represented by x, and the one at 500 K is The subtle 
difference in the middle lies in the range of thermal fluctuation 
at 500 K. 



On the other hand, the yield stress for an edge disloca- 
tion should be interpreted as the maximum stress during 
the whole pinning process including two partials: that is, 
the critical stress for the leading partial. It is interesting 
to compare the yield stress in the present simulation and 
the one in [l5l | in which the dissociation is not taken into 
account. The yield stress estimated by Scattergood and 
Bacon reads 
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which is also plotted in FIG. [S] We can see that Eq. © 
considerably overestimates the yield stress where r < 2.5 
nm. Since it does apply to bcc iron 0, H, 0, we can 
conclude that the dissociation makes depinning easier. 



C. temperature dependence 

In this section we investigate temperature dependence 
of the critical stress. We have calculated the critical 
stress for four different temperatures: 100 K, 200 K, 400 
K, and 500 K. However, we cannot find any differences 
regarding the critical stress between these calculations. 
In FIG.|n|we show the snapshots of dislocations just be- 
fore depinning where the temperatures are 100 K and 500 
K, respectively. We can see no difference in the disloca- 
tion shape. Also, the critical stress is almost the same: 
130 ± 12 MPa for r = 1 nm and 288 ± 30 MPa for r = 2.5 
nm. Namely, temperature plays no role in depinning pro- 
cesses. 

Note that this result is opposite to the simulation on 
bcc iron with a copper precipitate 0. There, definite 
temperature dependence was observed from K to 500 
K. The reason for the difference lies in the dissociation 
nature of dislocation. In bcc iron where a dislocation 
does not dissociate, a perfect dislocation absorbs vacan- 
cies from a void (or the precipitate surface) and under- 
goes climb motion. The climb motion is remarkable for 
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larger voids where temperature dependence of the criti- 
cal depinning stress is observed. Note that smaller voids 
show less temperature dependence and the climb motion 
is weak there. It implies that the climb motion is essen- 
tial to the temperature dependence. On the other hand, 
no climb motion is observed in the present simulation. It 
is known that the climb motion is difficult in fee crystals 
since dislocations dissociate. Therefore, no temperature 
dependence is observed in the present simulation on fee 
copper. 

Another possible reason lies in the activation energy 
for depinning. Although the precise estimation of the ac- 
tivation energy is difficult, it is at least larger than the 
energy of a dislocation whose length is equivalent to the 
void diameter. (We neglect the step formation energy on 
the void surface.) Dislocation energy is calculated by the 
line tension multiplied by the length. The effective line 
tension is estimated to be 0.4 nN by Eq. (|12|l . (Please see 
the next section for the detail.) Therefore, for the void 
of 2.0 nm radius, the activation energy is approximately 
estimated as 1.6 x 10~^^ J. It is equivalent to 400 fc^T, 
where denotes the Boltzmann constant and T ~ 300 
K. Since this is much larger than the thermal energy of 
involved atoms on the void surface (less than a hundred), 
it is plausible that thermal fluctuations cannot assist dis- 
location depinning. In addition, please recall that the 
step formation energy is neglected and the actual value 
may be larger than that. However, in the case of pre- 
cipitates, the number of involving atoms is larger than 
in voids. That may be another reason for the difference 
between the present simulation and the one of Q. 



IV. THE DISLOCATION SHAPE AND THE 
CRITICAL ANGLE 

A. an orientation dependent line tension model 

As can be seen in FIG. |31 the bowing dislocation is 
asymmetric with respect to a; = 0. Since the uniform 
line tension model predicts the symmetric form (an arc), 
we have to incorporate an orientation dependent line ten- 
sion in order to explain the asymmetric dislocation shape. 
Indeed, de Wit and Koehler obtained a solution for 
the dislocation shape. 
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where 5 denotes the angle between the tangent line of 
a dislocation and the x axis. An orientation dependent 
line tension is denoted by E{9), in which 9 represents 
the angle between the tangent line of the dislocation and 
the Burgers vector. If the concrete form of E{9) is given, 
these equations can be numerically solved with an appro- 
priate choice of the arbitrary constants Ci and C2. Here, 
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FIG. 7: Bowing dislocations obtained by the simulation for 
various voids (a: r = 1.0 nm, b: r = 1.5 nm, c: r = 2.0 nm, 
d: r = 2.5 nm.) The solid lines represent Eqs. O and JHJ, 
which are derived from an orientation dependent line tension 
model. Note that the trailing partials also deform to a certain 
degree as the void becomes larger. 



E{9) is given as [3 



(9) 



f{9) ~ 59.3- 16 cos 261 -0.8 cos 46*, [GPa] (10) 

where ro is the (arbitrary) core cutoff length. In FIG. [7| 
Eqs. © and H10|l are fitted with the simulation result. 
We can see that the fittings are quite satisfactory. 



B. determination of the critical angle 

In macroscopic materials, random configuration of ob- 
stacles plays an important role in the dislocation dynam- 
ics, as was discussed in the section I. In order to apply 
Eq. to a practical situation, we wish to estimate the 
critical angle from the simulation. Note that the critical 
angle 0c is defined by the angle between two tangent lines 
of the dislocation at the void surface. Because the critical 
angles are slightly different for two partial dislocations, 
we measure the both. 

The pinning strength a — cos(0c/2) with respect to 
the void radius is shown in FIG. |H1 We remark that the 
pinning strength a again obeys a logarithmic law. 



A log 



2r 



S(l + f)' 
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The constants are A — 0.24 and B — 0.07 for the leading 
partial, and A — 0.28 and B = 0.15 for the trailing par- 
tial. By extrapolation, a reaches 1 when the void radius 
exceeds 3 nm, as shown in FIG. |S1 Note that the ten- 
dency has also b_een observed in the previous simulations 
on bcc iron 
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FIG. 8: Void radius dependence of a for the leading partial 
(x) and for the trailing partial (+). 
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FIG. 9: Effective line tensions estimated via Eq. (inj. Those 
of the leading partial and of the trailing partial are repre- 
sented by X and +, respectively. Note that the void radius 
dependence results from the difference in the configuration of 
dislocations: the degree of bowing and the interactions of two 
partials. 



Meanwhile, it should be remarked that a crossover is 
not observed in the pinning strength a. From Eq. J^l, 
the crossover of should be attributed to that of 7, 
which can be written as 



(12) 



The calculated line tensions are shown in FIG.|2| in which 
the line tension of the leading partial shows the minimum 
at r = 1.0 nm while that of the trailing partial monoton- 
ically decreases. The explanation for the decrease lies 
in the deformations (i.e. bowouts) of dislocations. As 
the void gets larger from r = 0, the extent of bowout 
becomes greater. This results in the prevalence of the 
screw component, since the partial dislocation is initially 
edge-like. The deformation lowers the effective line ten- 
sion: i.e. the energy per unit length. (However, the total 
energy increases due to the elongation.) The above mech- 
anism also explains the decrease of the line tension of the 




FIG. 10: The impact parameter d is defined by the distance 
from the void center to the glide plane. Note that the lower 
part of the void corresponds to the negative values of d. 



leading partial for r < 1 nm. 

On the other hand, the increase of the line tension of 
the leading partial for r > 1 nm results from the in- 
teraction with the trailing partial. For larger voids, the 
trailing partial is also pinned by the void and bends dur- 
ing the depinning process of the leading partial, as can be 
seen in FIG.[7| That effectively increases the line tension 
of the leading partial. 



V. EFFECTS OF THE IMPACT PARAMETER 

So far, we have limited ourselves to the situation where 
a dislocation penetrates the void center. This is rather a 
special case, because the relative position of a void to a 
glide plane may be arbitrary. In this section, we change 
the distance between the void center and a glide plane. 
We call the distance "the impact parameter", which is 
denoted by d. (See FIG. [101) 

The pinning strength a is determined for the various 
impact parameters. The result is shown in FIG. 1111 which 
shows asymmetric dependence of a on d. Note that the 
dislocation is pinned even when it is not in contact with 
the void. It implies that, as well as the core energy, the 
elastic strain around the dislocation plays an important 
role in pinning processes. In addition, the asymmetry 
regarding d — comes from the nature of strain field 
around an edge dislocation: i.e. existence of hydrostatic 
pressure caused by the extra atomic plane. Especially, 
the fact that pinning strength for d > becomes con- 
siderably weak suggests that the hydrostatic pressure is 
dominant over the shear stress. 

Also, we remark that strong pinning (a > 0.5) occurs 
only where — l.Or < d < 0: i.e. the lower half of the void. 
This area accounts for approximately 30 or 40 percent 
of the whole pinning region, while the rest involves rel- 
atively weak pinning. The large variance of the pinning 
strength distribution for a single void suggests the recon- 
sideration of the same pinning strength assumption in 
dislocation dynamics simulations. In the section FVIII we 
discuss how to incorporate this effect into the estimation 
of the CRSS in the framework of Eq. 0. 
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FIG. 11: Impact parameter dependence of the pinning 
strength a for the traihng partial. Void radius is 1.0 nm 
(indicated by the dashed lines). 



VI. EFFECTS OF VOID DEFORMATION 
AFTER THE PASSAGE OF SEVERAL 
DISLOCATIONS 

When a void is sheared by a dislocation, two parts 
which are divided by the ghde plane are displaced to 
each other by the Burgers vector. After the passages of 
several dislocations, it may collapse and lose the pinning 
ability. For example, collapse of stacking fault tetrahedra 
by the passage of dislocations is both experimentally E3| 
and computationally ,5| observed. This phenomenon is 
believed to be responsible for the formation of disloca- 
tion channel and the localization of plastic flow, which 
recently invokes attention including some computational 
studies 0. In this section, effects of the void defor- 
mation on pinning strength are discussed based on the 
motivation described above. 

We remark that another possible mechanism of the 
void deformation is vacancy absorption by (and the climb 
motion of) an edge dislocation, as was discussed in the 
section Fill CI However, no climb motion was seen in our 
simulations, because the climb is difficult in fee metals 
due to its dissociation. Hence, we do not consider the 
void contraction by the vacancy emission to dislocations. 
We concentrate on the effect of the relative deformation 
with respect to a glide plane. 

We prepare the deformed void as shown in FIG. El 
First we prepare a spherical void. Then, instead of iter- 
ating the pinning simulations, atoms located above the 
glide plane are displaced by the Burgers vector a/2[110]. 
Iterating this procedure for TV times is equivalent to the 
passage of N edge dislocations. We set two configura- 
tions. In the case (a), the glide plane on which the dislo- 
cations move is assumed to be at the void center, whereas 
glide planes are uniformly distributed in the case (b). 

In the case (a), we set iV = 5 and N — 10. Namely, 
the void is assumed to be penetrated by dislocations on 
the same glide plane 5 times or 10 times. The pinning 
strength a is found to be 0.7 for the both cases, which 
is almost the same value with the one for the spherical 
void. This is consistent with the result obtained in the 



FIG. 12: Schematic of the void deformation by the passages 
of TV edge dislocations: a) There is a single glide plane which 
cuts the void center, b) Glide planes are uniformly dis- 
tributed. Each planes has dislocations which are to pene- 
trate the void. 



last section that the upper part of the void is dominant 
in pinning of edge dislocations because of the hydrostatic 
pressure. In the case (b), we test N = 2 and cannot find 
any difference from the non-deformed case regarding the 
critical angle. Thus, as far as the vacancy absorption 
mechanism (i.e. climb) is absent, the pinning strength is 
not seriously altered by the passage of dislocations. 



VII. DISCUSSIONS AND CONCLUDING 
REMARKS 

A. practical applications 

Let us estimate the yield stress of copper in which voids 
are randomly distributed. First, we determine the aver- 
age spacing L between obstacles on a glide plane. The 
areal density is represented as 2rp, where p denotes the 
number density of voids per unit volume. Then the av- 
erage spacing on a glide plane is written as L = l/y^2rp. 
Since the expression includes all voids which intersects a 
glide plane, their impact parameters are randomly dis- 
tributed from — r to r. Since we have seen that the pin- 
ning strength a considerably changes with the impact 
parameter d in the section^ we wish to incorporate this 
result. However, there is no simulation which considers 
this effect at this point. Therefore, we have to resort a 
rough approximation here. 

From FIG. ^2 strong pinning (a > 0.4) occurs only 
where — 1.2r < d < 0.2r. We neglect the rest. Namely, 
it is assumed that only this region is responsible for 
pinning. We use 1.4r instead of the diameter 2r; then 
L ~ 1/ y^lArp. We can rewrite Eq. |2l as 



1.87a 



(13) 



where the line tension 7 should be interpreted as the ef- 
fective one that is determined by the present simulation, 
FIG.H 
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For example, an irradiated copper specimen includes 
voids whose average diameter 2r = 4.1 nm, and the num- 
ber density p is 2.9 x 10^^ m"'^ |23|. From the present 
simulation, the line tension 7 and the pinning strength 
a are estimated as 0.42 nN and 0.9, respectively. The 
effective areal density is then calculated as 8.4 x 10^^ 
nm~^ (i.e. the average spacing is 110 nm). This yields 
Tc = 40 MPa, which is quite a reasonable value. In order 
to give more precise predictions, Eq. 10) should be mod- 
ified to include the effects of the impact parameter and 
the orientation dependent line tension. 



However, the problem is the arbitrary parameter B. 
This arbitrariness may vanish if we can calculate the 
deformed (curved) effect on the dislocation energy, al- 
though it is a very complicated problem. At least, we 
confirm that the extended form of the proportional con- 
stant Eq. (|17|l agrees with the present simulation, if there 
is no strong interaction between the partials. 



conclusion 



B. comparison with a continuous model with 
self-interaction 



Scattergood et al. have presented the calculation on 
the dislocation pinning by a void based on the frame- 
work of Bacon et al. jl3 • It is a continuous model that 
incorporates the self-interaction of a dislocation. Since 
their system also consists of the periodic array of voids, 
we wish to compare their result to ours. 

Let us briefly review their discussion. They have spec- 
ulated that depinning concerns the line tension near the 
pinning point, where a dislocation forms a dipole. Since 
a dislocation rotates by almost 7r/2 there, we regard the 
effective line tension as that of its counterpart. For exam- 
ple, as for an edge dislocation, the effective line tension 
upon depinning is that of a screw dislocation. Then the 
depinning stress for an edge dislocation is expressed by 



7off 
1 



27cff 
Lb 

47r 
B 
L-2r 



\ogR, 



B 
2^' 



(14) 
(15) 
(16) 



where R should be interpreted as the effective outer cutoff 
divided by the core cutoff. Note that B itself is not the 
core cutoff but an unknown function of the core cutoff. 
The above equations yield Eq. 

Their discussion can be extended to dislocations of ar- 
bitrary orientation. 



7off 



47r(l - 1^) 



1 2 / " 

1 - i^COS (- 



logi?, (17) 



where ly and 9 denote Poisson's ratio and the angle be- 
tween the dislocation line and the Burgers vector, respec- 
tively. In the present case, substituting 9 = tt/3 (for the 
partials) into Eq. (dTj) yields r = 30 log .R [MPa]. Recall 
that Eq. 1)17(1 is an expression for the resolved shear stress 
with respect to the Burgers vector of a partial disloca- 
tion. It is equivalent to t = 35 log .R [MPa] with respect 
to the Burgers vector of an original perfect dislocation. 
Indeed, it shows an excellent agreement with the critical 
stress of the trailing partial. 



We calculated the critical stress and the critical pin- 
ning angle for the interaction between an edge dislocation 
and a void in fee copper. Dissociation of a dislocation 
plays an important role in the behaviors of the critical 
stress: i) It is much lower than the estimation of Scatter- 
good and Bacon, which does not consider dissociation, ii) 
It suddenly increases at a certain void radius where two 
partials can be simultaneously trapped, iii) The depin- 
ning stress of the trailing partial does agree with that of 
Scattergood and Bacon, since the leading partial moves 
far from the void. 

We also found that there is no temperature dependence 
in the critical stress and the pinning angle. This is op- 
posite to the previous simulation on bcc iron 0, H, 0- 
The difference comes from the presence (in bcc) or the 
absence (in fee) of climb motion. 

The pinning strength cos(0c/2) obeys the empirical 
logarithmic law which has been found in 0, . The 
distance between the void center and the glide plane (the 
impact parameter) is found to affect the pinning strength 
in asymmetric manner. This is due to the hydrostatic 
pressure around an edge dislocation. Hence, it is inter- 
esting to compare the result with that of a screw dislo- 
cation, which is a work in progress. 

The impact parameter dependence of the critical an- 
gle also suggests the importance of randomness in the 
pinning strength. Even if the system contains voids of 
the same radius, the cross section on a glide plane is 
randomly distributed. Hence we have to incorporate the 
randomness in the pinning strength. It is not straightfor- 
ward to deduce this effect from the existing simulations 
which treats only two kinds of obstacles [22]. Investiga- 
tion of a continuous model with random pinning angles 
will be interesting to see how the impact parameter de- 
pendence affects the macroscopic dislocation motion. 
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